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Abstract 

The problem of diagonalization of the quantum mechanical Hamiltonian, governing 
dynamics of an electron on a two-dimensional triangular or square lattice in external uni- 
form magnetic field, applied perpendicularly to the lattice plane, the flux through lattice 
cell, divided by the elementary quantum flux, being rational number, is reduced to gen- 
eralized Bethe ansatz like equations on high genus algebraic curve. Our formulae for the 
trigonometric case, where genus of the curve vanishes, contain as a particular case recent 
result of Wiegmann and Zabrodin. 
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1. Introduction 

In this paper we consider the diagonahzation problem of the following Hamilton operator: 

n = iJi{aS + a-^S-^) + v{(3T + (i'^T-^) + p{-iU + 7"'f/"'), (1-1) 

where unitary operators T,and U satisfy Weyl commutation relations 

ST^iuTS, TU = iuUT, US = ujSU (1.2) 

with u being a primitive A'"-th root of unity: 

a; = exp(27riM/iV), (M, N) = 1 (1.3) 

for some mutually prime integers AT > M > 1, complex parameters a,/?, 7 have unit 
absolute values, while fx, v and p are real parameters. Phases of a, 7 are constrained to 
lay betweem and 2t^/N: 

< arg(a), arg(/?), arg(7) < 2tt/N. (1.4) 

Hamiltonian (1.1) appears in various problems of solid state physics. With parameter p, 
taken to be zero, it governs the dynamics of an electron on a two-dimensional square lattice 
in the external uniform magnetic field, applied perpendicularly to the lattice plane, phases 
of a and f3 playing the role of quasi- momentum, and the phase of a;, the role of magnetic 
flux through lattice cell, divided by the elementary quantum flux. The same problem on 
triangular lattice is described by (1.1) with non-zero parameter p. Physical meanings of 
parameters //, v and p are length scales along corresponding directions on the lattice. 
A'"-th powers of operators S, T, U are central elements and we flx them to be unity: 

5^ = = t/^ = 1. (1.5) 
Prom (1.2) and (1.5) one can deduce that operator 

C = STU (1.6) 
is also central element, and its A'"-th power is 1 (—1) for odd (even) N: 

C^ = (-1)^-1. (1.7) 

If we fix C by some number, satisfying (1.7), then operators S,T,U have a unique A'"- 
dimensional irreducible representation. Por the problem under consideration only algebraic 
properties of these operators, given by (1.2), (1.5) are important. Different reducible 
realizations can affect degeneracies of eigenvalues, but not the spectrum itself. We will use 
this fact, working with A^*^ -dimensional reducible representation. 

Exploration of the spectrum of (1.1), or its various particular cases, is associated with 
many names, and especially Hofstadter, who at p = made a detailed numerical analysis 
of dependence of this spectrum on integers M and N ior N < 50 [H] . Shortened list 
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of physical and mathematical papers, concerning this problem, contains also [W], [HK], 
[BKS], [CEY]. 

Last years the idea about possible relevance to above problem of methods, used in 
quantization of integrable models, attracted some attention. This idea has been mate- 
rialized in recent paper [WZ], where Wiegmann and Zabrodin showed that in the case, 
equivalent to p = 0, /i = —v = 1, and = (3^ = 1, the spectrum of (1.1) can be 
represented in terms of solutions of Bethe ansatz (BA) type algebraic equations. 

In this paper, using more sophisticated methods of integrable models, we show that 
BA like equations can be obtained also for general Hamiltonian (1.1). They are similar 
to the equations, found in [BS] for the case of the chiral Potts model. Rapidity variables, 
satisfying generalized BA equations, live on high genus algebraic curve. The rational limit, 
where the genus vanishes, contains also the case, considered in [WZ]. We hope, that our 
result will be useful for further investigation of Hamiltonian (1.1), especially in the limit 
of large A^. 

In Section 2 we represent Hamiltonian (1.1) as a part of three-site transfer matrix, 
constructed through elementary L-operator, intertwined by the six-vertex i?-matrix. In 
Section 3 we calculate Baxter's vector, which enables to derive functional equation, de- 
termining in principle the spectrum of Hamiltonian (1.1) in terms of solutions of BA like 
equations on high genus algebraic curve. Section 4 contains specialization to the case with 
p = in (1.1). In Section 5 a detailed investigation of the trigonometric case, corresponding 
to genus zero curve, is performed. 

2. Transfer Matrix 

In this section we introduce a simple L-operator, intertwined by the six-vertex i?-matrix, 
and define transfer matrix, containing Hamiltonian (1.1). 

Consider two operators, X and Y, satisfying Weyl permutation relation, and with 
unit A'"-th powers: 

XY = u;YX, X^ = Y^ = 1, (2.1) 

with oj from (1.3). These operators can be realized as N-hy-N matrices with the following 
matrix elements: 

{m\X\n) ^ u'^dm,n, {mlYln) ^ dm,n+i, (2.2) 

where indices m and n run over N values 0, 1, . . . , A" — 1 and are considered (mod N), 
so they are elements of Z^. 

Introduce 2-by-2 matrix, L{x, h), with operator valued matrix elements: 

^(-'•)=("cl T)' ^'-'^ 

where x and h = (a, b, c, d) are some complex parameters. The A/"-dimensional linear 
space, where operators X and Y act, will be refered to as a "quantum" space while the 
two-dimensional space, where L(x, h) acts as 2-by-2 matrix, as an "auxiliary" one. 
Matrix (2.3), called L-operator, is intertwined by the six- vertex i?-matrix: 

R{x/y)L{x, h) ® L{y, /i) = (1 ® L(y, h)){L{x, h) ® l)R{x/y), (2.4) 
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where L{x, h) and L(j/, h) act independently in two different auxiliary spaces and in one and 
the same quantum space, while R{x/y) is a matrix in the tensor product of the auxiliary 
spaces with numerical matrix elements: 



R{x) 



(xu-x-^ \ 

uj{x-x-^) a;-l 

a;-l x-x-^ 

V xoj-x-^) 



(2.5) 



Note, that this -R-matrix differs from the usual one in two diagonal elements. Our choice 
enables us to escape square roots of uj for the time being and to consider both odd and even 
N simultaneously. Matrix (2.3) has been used in [BBP] for derivation of functional equa- 
tions for the chiral Potts model. It can be extracted also by specialization of parameters 
in L-operators, written in [BKMS] and [T] . 

Now introduce the following transfer matrix, acting in the tensor product of three 
quantum spaces: 

T [x) = tr(L(x, ho) O L{x, hi) ® L{x, /is)), (2.6) 

where the matrix products and trace are performed in the auxiliary space, and the indexed 
h's mean that three different sets of parameters hi — {ai, bi,Ci,di) are taken. To write down 
explicitly T{x), introduce special notations for certain combinations of parameters hi: 

Ci = bi-idiCi+i, fi = Ci-ittibi+i, z = 0, 1,2 (mod 3), (2.7) 

where index i is considered as an element of Z3. Then, for T{x) we have the following 
explicit expression: 

r{x) = To + x''T2, (2.8) 



where 
and 
with 



% = ojaoaia2C + dodid2, (2.9) 
% = €28 + + eiT + /iCT-i ^ ^^jj ^ f^cU-\ (2.10) 

S = X®Y®1, T = y®l®X, U = 10X0Y, (2.11) 

C = uj-^XY (g)XY (g)XY. (2.12) 

Operators S,T,U and C, defined in (2.11) and (2.12), satisfy algebraic relations (1.2), 
(1.5), and (1.6), operator C being commutative with others (so it can be considered as a 
number). Thus, operators (1.1) and (2.10) can be identified: 

T2 = n, (2.13a) 

where 

a = el/\f2C)-y', ^ = e'l\hC)-'l\ l = e]!\m-'l\ i^lZh) 
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/^=(e2/2C)V2, ^ = (e^j^c')i/2^ (2.13c) 

To conclude the section note, that one and the same transfer matrix T{x) can be 
represented by (2.6) with different L's: 

r (x) = tr(Lo(a;, ho) Li{x, hi) L2{x, /ta)), (2.14) 

where L's differ from L's by a gauge transformation: 

Li{x, hi) = AiL{x, hi)A-^^, i e (2.15) 

with some invertible numerical 2-by-2 matrices Ai. In the next section we will use this 
freedom to construct Baxter's vector. 

3. Baxter's Vector 

Following [B] , [FT] , and [BS] , let us turn to the calculation of Baxter's vector for transfer 
matrix (2.6). For this let us use its representation by (2.14) with matrices Ai, chosen as 

I l), Vi = ^i-h i&Zs, (3.1) 

where are some complex parameters to be fixed later on. Transformed L-operators 
Li{x, hi) have the following explicit form: 

f /, ^ _ -F{x,hf,r]i,r]i+i)\ , . 

^'^""^^'^ - \F{xMMM-,i) -F{x,hi;^i,r}i+i))^ ' ^ ^^'^^ 

where operator valued function F is defined by 

F{x, h; e, e') = i^'aX - xb)Y + ^(x^'cX - d). (3.3) 

Now let us try to fix parameters by demanding, that equations 

F{x,hi;$i,^i+i)\(l>i) = 0, ieZs (3.4) 

have non-zero solutions. Writing them in the basis (2.2), we get the following recurrence 
relations for matrix elements of vectors 10,): 

= -^i It '..m ^eZs. (3.5) 



(m - l|0i) * x^i+iCjo;'^ - d. 
Periodicity conditions impose algebraic relations of high order on ^'s: 

- (-1) cf (..Cm)^ - cif ' ^ ^'-'^ 
4 



Considering ^'s and x as variables, while keeping all /I's fixed, we get three equations on four 

variables, which define an algebraic curve F. Let p be some point of this curve, specified 
by a particular set of coordinates (x,^^), satisfying (3.6). Then vectors depend on this 
point, so from now on we will use another notation for them: 

\(l>i) = \p)i, per. (3.7a) 

Besides, choose the following normalization: 

{0\p)i = 1, ie Zs. (3.76) 

To proceed further, define two automorphisms of our curve, t±, which act on coordinates 
X and as follows: 

T±:r^r, rix = a;±i/2x, r^^ = ^"'/'ei- (3-8) 
Using this definition, one can easily show that 

F{x,hi;^i- l,^i+i)\p)i = -\T_p)i{x^i+iCi- di), (3.9a) 

and 

F{x,hi;^i,^i+i - l)\p)i = -\T+p)i^i{aidi - x%iCi) / {^i+^ai - xbi). (3.96) 

Denoting 

b) = b)o<H)b)i®|p)2, (3.10) 

and using (2.14), (3.2), (3.4), (3.9) together with multiplication properties of triangular 
matrices, we come to the relation: 

T{x)\p) = \t_p)A_{p) + \t+p)A+{p), (3.11) 

where 

A_{p) =Y[{di- x^i+^a), (3.12a) 

^+{P) = n ^ii^i^i ~ xHiCi) / {^i+itti - xbi). (3.126) 

To get the functional equation from (3.11) we just multiply it from the left by the 
eigenvector {(p\ of T{x), corresponding to some eigenvalue A(a;). Vector {(p\ evidently does 
not depend on p, so we obtain eventually the scalar relation 

A{x)Q{p) = Q(t_p)A_(p) + Q{t+p)A+{p), (3.13) 

where 

Q{p) = {<P\P) (3.14) 
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is a function on our algebraic curve with known poles (they can be extracted from recur- 
rence relations (3.5)), while its zeros, the number of them being equal to that of poles, are 
determined from the generalized Bethe ansatz equations, obtained from (3.13) by taking 
p as various zeros of Q{p): 

/, = l,...,#(poles). (3.15) 

Q[r+Pk) ^-[Pk) 

Relations (3.13) - (3.15) at this stage only in principle solve the diagonalization problem of 
(1.1), since eqs. (3.15) are too complicated to work with. The main problem, of course, is 
that of the suitable choice of coordinates on F. Unfortunately, uniformization of algebraic 
curves is still unsolved in general, long standing problem of mathematics. Our hope, 
however, is that F has a very special structure, which we believe should admit particular 
approach to it. Anyway, we think that the analytic structure of F deserves a further study. 

4. Hofstadter Hamiltonian 

Here we specify results of Section 3 to the case with p = in (1.1), which corresponds to 

ao = do = 0, 6o = Co = 1- (4.1) 

The algebraic curve F, defined by (3.6), reduces to a disjoint set of N copies of another 
curve, Fo, defined by two high order equations on three variables ^0)^2 and x: 

NtN _ uN N NtN _ hN N 

^° cf(xe2)^-<' ^ c^ix^or-d^' ^^-'^ 

while s-th copy in the set is specified by the value of ^i: 

Ci = ^'-'/'/eo, s e Zn (4.3) 

with — o;^/^ being chosen as A'^-th root of minus unity: 

Now Baxter's vector \p) as well as functions A±{p) acquire index s: 

\p) \p, s) , A± (p) A± (p, s) , seZn (4.5) 

so relation (3.11) now reads 

T{x)\p, s) = |t_p, s - l)A_(p, s) + \t+p, s - l)A+{p, s), (4.6) 

where automorphisms t± of Fq are defined by (3.8) with index i being restricted to only 
two values, and 2. Taking into account (4.3) together with (3.12), we have explicitly: 

A±(p,s)=a;«A±(p,0). (4.7) 
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Noting that 

a;" = $(s)$(l)/$(s-l), = a;*(^+^)/^ (4.8) 

then multiplying (4.6) by a;**/$(s) and summing over s, we get Baxter's relation for the 
case of Hofstadter Hamiltonian: 

r(x)|p,t)'= |T_p,t)'A_(p,t)$(l) + |T+p,t)'A+(p,t)$(l), teZM (4.9) 

where 

\p,t)'= J2 \p,s)u''ms), (4.10) 

sGZn 

and we used (4.7). Repeating arguments from the end of Section 3, one can get direct 
counterparts of relations (3.13) - (3.15). 

5. Rational Limit 

In this section we consider rational limit of formulae from Section 3. 
Let us restrict parameters hi, i & Z^, by 

ai = q~^di, hi = q~^Ci, i e Z3. (5.1) 

Everywhere in this section q = u?-!"^ satisfies (4.4). Then, relations (3.6) become very 
simple: 

^f = e^ ^eZ3, e'^ = l. (5.2) 

Here variable x does not enter, so, we have just several copies of genus zero curve, spanned 
by coordinate x, while all ^'s are fixed up to A'"-th roots of unity. Let us put ^ = q^ for 
some 1 = 0,..., 2N — 1 and choose 

= i^Zs, (5.3) 
then, point p of F can be identified with pair (x, I), and (3.11) becomes 

T{x)\x, I) = \xq~^, I - l)A_(a;, I) + \xq, I - l)A+(a;, I), I = 0, . . . ,2N - 1, (5.4) 
where 

A_(x,l) = Ylidi- Cixq^), (5.5a) 
A+{x, I) =11 {d1 - clx'')/{di - axq-'). (5.56) 

Define function 

[1/2] 

fix, l)= JJ n - xCiq-^Lo')/{d^ - xcq^u-'), I = 0, . . . ,2N - 1, (5.6a) 



where 



. / ._n/2, / = (mod 2); 

(/-l)/2, 1 = 1 (mod 2), ^^-^^^ 



and for any m & Zn introduce new vectors 

N-l 



x,m)e= ^2 \x,2n)f{x,2n)<jj^^, 

n=0 
N-l 

X, m)o = \x, 2n + l)f{x, 2n + l)a;'^ 



(5.7) 



n=0 

Combining them into 2-component row vectors 

\x,m)) = {\x,m)e \x,m)o), (5.8) 
one can rewrite (5.4) in a matrix form: 

T{x)\x,m)) = \xq~^ ,m))D-{x,m) + \xq,m))Dj^{x,m), (5-9) 



where 



Relations (5.9) do not change their form under gauge transformations: 

\x,m)) \x,m))U{x,m), 
D±{x,'m) — > U{xq^^ ,'m)~^ D±{x,m)U{x,'m) 



(5.11) 



for any invertible 2-by-2 matrix U{x,m). In the case of arbitrary and without fur- 
ther restrictions on parameters, there is no a gauge, where matrices D±{x,m) would be 
diagonal, and rational in x, so, from now on let us work with only odd N: 

N = 2P+1, P>1. (5.12) 

In this case q = w^/^, satisfying (4.4), can be represented as an integer power of to: 

q = ^i/2=^P+\ (5.13) 

Define function 

p 

u{x) = II l[{di - xciu;^), (5.14) 

ieZa j=0 

and choose gauge transformation matrix in (5.11) of the form 

TTt \ f uiqx) q~'^u{x)\ 



Then matrices D± {x, m) become diagonal: 

i;±(x,m) = rV'(±^9^'/')(j (5.16) 

where 

il}{x)= Y[{di + xciq-^^''). (5.17) 

Now multiplying (5.9) from the left by the eigenvector {(f \ of T(x)^ corresponding to some 
eigenvalue A(a;), we get two scalar relations: 

±q-'^K{x)Q±{x) = ^(-xg-i/2)Q±(xQ-i) + il^{xq^'^)Q±{xq), (5.18) 

where 

(g+(x) g_(x)) = (^|x,m)). (5.19) 

From the structure of Baxter's vector, given by (3.5), it follows that Q±{x) is a polynomial 
in X, and sending x to zero and infinity, taking into account (2.8) and (2.9), we conclude, 
that 

Q_(a;)=0, Q+{x)^x^Q{x), degQ{x) = 2P, m = P, (5.20) 

while central element C in (1.6) equals to unity. For convenience choose the following 
normalization for the parameters cq, ci, C2*: 

C0C1C2 = w- (5-21) 

Then, function ip{x) reads 

ilj{x) = q^/^{x + n){x + iy){x + p), (5.22) 

where jJ^ji^jP are the same as those in (1.1). As for the parameters a,/?, 7, they are equal 
to q^^'^. Using (5.20), we rewrite (5.18) as 

A{x)Q{x) = i;{-xq~^/^)Q{xq~^) + q~^i;{xq^/^)Q{xq), degQ{x) = 2P, (5.23) 

A (a;) being 

A(x) = iiupiq^/'^ + q-'^/'^) + x'^E, (5.24) 
where E is an eigenvalue of Hamiltonian (1.1). Let Q{x) has the following decomposition: 

2P 

Q{x) = n (^ - V-^m) (5.25) 

m=l 



* This can be done by rescaling variable x 
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for some z^s, which should satisfy BA equations (3.15): 



-1/2 {^izi + q^/^){uzi + q^/^){pzi + q^/^) |^ qzi - Zm . . 

^ {q^^fizi - l){qy^uzi - l){q^/^pzi - 1) }} ^zi - qz^' ^"^^ ^""'^^ 

Differentiating (5.23) twice with respect to a; at a; = 0, we get the expression for E entirely 
in terms of z^s: 

E = fxMq-q-^){q^/^-q-^/^) ^ ZmZn 

l<m<n<2P 

2P (5.27) 

-(g - Q~^){i^^ + ^p + PiA 5^ + (g'/' + g"^/^)(/i + p + p). 

m=l 

Formulae (5.26) and (5.27) at p = 0, /U = — z/ = 1 reproduce correspondingly formulae (6) 
and (5) of paper [WZ], if we identify our Zm with their izm- ** 

It is worth to note that in the case of odd N one can start from the very beginning 
with another L-operator: 

where operators V and W satisfy 

VW = qWV, = W^ = 1. (5.29) 



L-operators (5.28) and (2.3) are connected in a simple way. If we multiply (5.28) from the 
right by V, and replace parameter a by qa, then we get (2.3) with 

X = W-^V, Y = WV. (5.30) 



L-operator (5.28) can be extracted as a particular case from the L-operator, considered 
by Bazhanov and Stroganov in [BS]. The case, where a = d and b = — c, is related to the 
massless quantum Sine-Gordon model on the lattice [G], [V], [KT]. 

Using L-operator (5.28), we can repeat the procedure developed in Sections 2 and 3. 
But automorphisms t± in this case act on variables and x in another way, namely, they do 
not change ^'s at all. So, in the rational limit we are not forced to use simultaneously several 
disjoint copies of the zero genus curve. Particularly, (5.23) and (5.26) arise as a direct 
specialization of general formulae, counterparts of (3.13) and (3.15), without intermediate 
transformations, which we made in this section. Nevertheless, the approach which we have 
presented in detail, enables us to consider both even and odd N on equal footing. This 
explains why we prefer to work with L-operator (2.3). 



In fact our eqs. (5.26) differ from (6) of [WZ] by a sign. We checked our equations by 
direct diagonalization at A/" = 3, so we think that there is a misprint in [WZ]. 
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Summary 

The main result of this paper consists in estabhshing the relation between the generalized 
Hofstadter Hamiltonian (1.1) and inhomogeneous XX Z chain of three sites. The available 
methods in lattice integrablc models lead to the generalized Baxter relation (3.11), as well 
as Bethe ansatz (BA) type equations (3.15) on high genus algebraic curve. In particu- 
lar rational case we have got BA like equations (5.26) together with formula (5.27) for 
eigenvalues of the Hamiltonian (1.1). 

To use Baxter's relation in general case, one has to develop analysis on the high genus 
algebraic curve. Apparently, one needs more information about automorphic functions, 
connected with the curve. 

The tantalizing opportunity also is to try to apply the known simplification properties 
of BA equations when the number of variables, satisfying these equations, tends to infinity. 
In our case this corresponds to N ^ oo. However, at the moment it is not clear, how to 
utilize this possibility. 
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